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A MODEL FOR CHEMICALLY REACTING NITROGEN-OXYGEN MIXTURES 


WITH APPLICATION TO NONEQUILIBRIUM AIR FLOW 

By Walter A. Reinhardt and Barrett S. Baldwin, Jr. 
Ames Research Center 


SUMMARY 


An air model is presented that is relatively simple yet quantitatively 
realistic for use in computing chemical nonequilibrium in a flow field. Such 
a model is applicable to those computational problems associated with high- 
speed reentry where the effect of the nonequilibrium processes on the primary 
flow variables would be important, but where the details of these processes 
would be of secondary interest. The constituents that comprise the system are 
N 2 , 0 2 , N, 0, N + , 0+, and e“. The dissociation-recombination reactions of N 2 , 
0 2 , and NO as well as the bimolecular exchange reactions are considered. How- 
ever, the presence of NO is accounted for in an approximate manner to avoid 
coupling in the atom conservation equations. The model allows for the inclu- 
sion of either nonequilibrium or equilibrium vibration while ionization is 
taken to be in equilibrium. The resulting rate equations are in a form that 
may be integrated over local values of the flow variables to obtain a system 
of transcendental equations that can be solved by iteration. Comparisons with 
experimental nozzle flow results and with other numerical calculations are 
presented. 


INTRODUCTION 


Developments in nonequilibrium flow theory over the past 10 years indi- 
cate a need for two types of gas models. Air models based on simplifying 
approximations have been an important source of qualitative information, 
whereas attempts to describe specific experimental results have led to great 
complexity in the models considered. This report presents a simplified model 
for nitrogen-oxygen mixtures that contains only two nonequilibrium variables 
for the chemical effects and yet is quantitatively realistic over a wide range 
of flow conditions • 

Early treatments of chemical nonequilibrium in the flow of pure diatomic 
gases utilized the Lighthill-Freeman model (refs. 1 to 3 ). This method is based 
on the assumption that in the presence of chemical reaction any concurrent 
effect of vibrational nonequilibrium will be minor; a ,T half excited" value of 
vibrational excitation is used. Variations on the method allow for a more 
accurate evaluation of the electronic and vibrational excitation in the equi- 
librium constant and enthalpy (refs. 4 to 6 ) . For nitrogen- oxygen mixtures, 
it was determined on the basis of available rate constants that reactions 
involving nitric oxide would play an important role in shock tube flows 
(ref. 7 )* A similar reaction model including vibrational nonequilibrium was 


incorporated in a machine program for channel flow (refs. 8 and 9). We shall 
have occasion to use the results of Emanuel and Vincent! (ref. 9) for purposes 
of comparison in the present paper. An excellent summary of available infor- 
mation on vibrational and chemical reaction rates believed to be of importance 
in air flows at temperatures up to 8000° K was given by Wray (ref. 10) . A 
pertinent set of calculations on channel flows and external flows by 
Eschenroeder et al. , emphasize the role of the bimolecular exchange reactions 
involving nitric oxide (ref. 11 ) . The reaction model used was essentially 
the same as that of the previous papers except the vibrational excitation was 
taken to be in equilibrium (we shall also use a representative sample of these 
calculations as a comparison for our model results) . The implications of the 
existing knowledge on nonequilibrium effects in afterbody flows have been dis- 
cussed (ref. 12) and possibilities for simplifying the reaction model in 
certain regimes in studies of complicated geometric configurations were 
considered (ref. 11 ) . 

In the foregoing works a mechanism for coupling between vibrational and 
chemical nonequilibrium postulated by Hammerling, Teare, and Kivel (ref. 13) 
was suppressed in the belief that it would be unimportant except at higher 
temperatures . The available analyses of shock tube data supported this view. 
Subsequently the theory of such coupling was improved by Heims (ref. 14) and 
Treanor and Marrone (refs. 15 and l6) and used to explain the results of shock 
tube experiments in an argon-oxygen mixture at temperatures up to 16,500° K. 
Additional coupling mechanisms of possible importance suggested by Bauer and 
Tsang (refs. 17 and 18) include a rotation-vibration coupling effect and a 
kinetic effect, due to the bimolecular exchange reactions, that would tend to 
push the vibrational excitation more rapidly toward equilibrium. Two reports 
(refs. 19 and 20) have been published that contain a discussion of much of 
the recent work in the nonequilibrium field as well as details on a general 
machine program that was used in many of the foregoing studies . 

Recently evidence has been uncovered to support the view that (at least 
in expanding flows) the rates of vibrational relaxation are much greater than 
the values obtained from previous analyses of shock tube data (refs. 21 to 23). 
In air mixtures there is evidence that the bimolecular exchange reactions may 
also promote vibrational relaxation (refs. 17, 18, and 22). Consequently, the 
accuracy of nonequilibrium theory, as it applies to vibrational relaxation, is 
still subject to doubt pending further experimental and theoretical investiga- 
tion. Until the remaining uncertainties are resolved, it is considered worth- 
while to allow for the possibility of uncoupled vibrational nonequilibrium. 

In the following analysis it is a simple matter to include this option and 
this has been done. 


NOMENCLATURE 1 


a ratio of total mass of oxygen to the total mass of nitrogen in a 

sample of the gas mixture (0.3064 for air) 

b ratio of gram molecular weights, Mo/M^ (1.142) 

^gs. units are used. 
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e eZ 

®eZ,i 

e vi 

f N,fO 


g i,j 

h 

ft 

% 

k 

Mi 

m± 

P 

P r 

R 

R 

T 

T r 


internal energy contribution due to electronic excitation, equa- 
tion (5) 

electronic internal energy contribution due to species i, equa- 
tion (6) 

vibrational contribution to internal energy due to species i, 
equation (l8) 

variables, defined by equations (27) and (28), that are in the 
rate equations for production of nitrogen and oxygen 

degeneracy of the jth term level for species i 

specific enthalpy, equation (3) 

Planck's constant divided by 2 rt 

equilibrium constant for ith reaction 

Boltzmann 1 s constant 

forward rate constant of the ith chemical reaction 
gram molecular weight of species i 
atomic weight of species i 
pressure 


reservoir pressure 

electronic partition function for species 

gas constant, equation (7) 

universal gas constant 

temperature 

reservoir temperature 


i, £ Sij exp (-^i) 
J 


W N ,WQ 


z 


degrees of dissociation for nitrogen and oxygen, respectively 
degrees of ionization for nitrogen and oxygen, respectively 
compressibility, equation (2) 


P I,N* P I, 0 J 


variables, defined by equations ( 12 ) to (15), and contained in the 
equations for equilibrium amounts of w^ , wq , y-^ , and Yq 9 
respectively 
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G i 


£ 


ij 


specific mole concentration, moles of i per unit mass of gas 

efficiency factor of the Zth species for promoting the ith 
reaction 

measure of the departure of species i from local equilibrium, 
v ie “ w i 

electronic energy corresponding to the jth spectral term for 
species i 




characteristic temperatures for rotation, vibration, dissociation, 
and ionization, respectively 


P 


density 


p Q standard density 

t characteristic relaxation time 


T 


v i 


characteristic relaxation time for vibration associated with 
molecular species i, sec 


DESCRIPTION OF THE MODEL 


In a well-known treatise on the thermodynamic and transport properties of 
high-temperature air, C. F. Hansen (ref, 24) devised a simplified model for 
air in equilibrium. The part of this work dealing with thermodynamic prop- 
erties has fallen into disuse because of the feasibility of tabulating (and 
including in machine programs) the results of more accurate calculations. In 
the case of nonequilibrium flows, however, such tabulations are not practical. 
Consequently, the basic ideas of Hansen are useful in a nonequilibrium model 
for air and have been applied with slight modifications in the present work. 


Equations of State 

The species to be included in the expressions representing the equations 
of state are 

N 2 , 0 2 , N, 0, N + , 0 , e" 

At equilibrium any other species that may be present in trace amounts do not 
affect the equations of state importantly. We assume that this will be true 
in the nonequilibrium states of main interest also. The omission of NO from 
the system permits the use of Lighthill variables, that is, degrees of disso- 
ciation and ionization. The atom conservation equations are thereby removed 
from the system of equations to be solved in nonequilibrium problems. Then 
the equations of state can be written 


1 + 



p = RpZT 


( 1 ) 


1 + 


Z = 


+ 2y w + | (w Q + 2y 0 ) 


1 + ,- 


( 2 ) 


h = - RT + RT 1 + 
2 


3v n + 10y N + | (3w 0 + 10y 0 ) 


Dft[ 


I 


+ e 


el 


(3) 


where 
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(fe) 
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y N e el,N H 





y O e e l,0 + J 



'eZ ,i 


R 



2(1 + a)Mjj 



g i 
1 ><. 


€ ij 

kT 


exp 



(5) 

(6) 


(7) 


The term e v is the vibrational contribution to the energy (per unit mass) 
and is given by the sum of equations (4a) and (4b) when vibrations are in 
equilibrium. The description of e v for nonequilibrium molecular vibration 
is given subsequently. The parameter a is the ratio of the total mass of 
oxygen to the total mass of nitrogen in a sample of the gas mixture, and b 
is the ratio of molecular weights Mq/j% . The quantities w N and y N are the 
Lighthill variables representing the degrees of dissociation and ionization, 
respectively, of nitrogen; Wq and yQ are the corresponding variables for 
oxygen. The relationship between these variables and the species concentra- 
tions is given in the appendix. Table I contains the values of the fixed 
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constants. These expressions correspond to Hansen's model except for the use 
of truncated harmonic oscillators for the vibrational energies and the addi- 
tion of a few more terms in the electronic energies of the species. The local 
equilibrium values of the Lighthill variables can be expressed in the 
following form that is convenient for calculation. 


W N P 


v 0 = 


y^ + 

yo + b, o 


1 " y NV /2 

1+2 — — - - 1 


%),N 


^,0 


1/2 


1 - yo, 

1+2 — -1 - 1 


y Np “ 2 V Pl ^ + F y 0 J 


1 + 




N 


p i,n + 1 y 0( 


1/2 


i 


(8) 

(9) 

( 10 ) 


\ “II (|3 i,o + y Me ) 


1 + 


^ £ v QPl,0 

(Pi,o + yN e ) J 


1/2 


(ii) 


where 


- mivr0. 


1\ 3/2 -2. 

%i,x 2 


1 - exp 


-a 


v,N. 


T 


- _ b 1 + a - / m Q k \ 3 ^ 2 Qei,0 rz 

l%,0 " a 2p~ B R -°aWV a . . 


1 - exp 


1 - exp 


-0- 


D,0; 


exp 




T 


( 12 ) 


'el , 0 2 


-a 


v,0 : 


T 


1 - exp 


• 0 D,O- 


exp 


%o 2 


Pl,N = 


2(1 + a) /m e kT'\ 3/2 Se‘Z,N + 


m w 


N \2jih 2 / Q 


exp 






(13) 

(14) 

Pl,0 = ' p ' m N ^2 ) (15) 

and the Q , . are the electronic partition functions for species j . Again 
numerical values of the constants are given in table I. 


2(1 + a) /m e kT\ 3/2 %l, 0+ ^- 0 I,O‘ 


We have compared the values of compressibility Z and specific enthalpy 
h at equilibrium with those of Hilsenrath and Beckett (ref. 25) and found 
them to deviate by less than 3 percent at temperatures up to 25>000° K in the 
density range 0.001-1.0 times atmospheric density. 
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Ionization has been included in the equations of state so that calcula- 
tions involving high-temperature equilibrium regions can be made without 
switching to a different gas model# Ionization nonequilibrium will not be 
considered# Thus wherever the quantities y^ and y 0 appear (eqs. (2) to (5), 
( 8 ) and (9))^ they should be set equal to y^ and y 0 ^, respectively (eqs. (lO) 

and (ll)). At temperatures below about 10 , 000 ° K, y^ and Yq could be set 
equal to zero . 


Rate Equations 

Rate equations leading to values of e v , w^ and Wq appearing in equa- 
tions (2) and (3) are needed to complete the model. Ionization and electronic 
excitation are taken to be in equilibrium. 


Vibrational nonequilibrium equations .- The rate equations required for 
uncoupled vibrational nonequilibrium are in a form 

= -(tvi) ^vi + Av ± ( 16 ) 

where 


e vi - e vei “ e v± 


(17a) 




devei 

dt 


(17b) 


The subscript v denotes vibration, subscript e refers to local equilibrium 
value (eqs. ( 4 a) and ( 4 b)) and the i indicates species (N 2 or 0 2 ) . The 
variable t y ± is a characteristic relaxation time (values of t y ± obtained 
from ref. 19 are given in table 11 (b)). The quantity e v in equation (3) is 
to be evaluated as 

e v = e V]fl 2 + e VQ 2 ( 1 2 3 ^) 


rather than e y = e v e 0 2 + e VeN 2 as i n( ii ca “ked previously for vibrational equi- 


librium. 

example. 


The quantities 


e V]y 2 aad e v Q 


e Vi\- 


2 

= °v e N 2 


are obtained from equation (17a); for 

(19) 


V W 2 ~ -v N2 

Chemical nonequilibrium equations .- The chemical equations of main inter- 
est are the dissociation reactions 


1 . K 2 + M ^ 1 + N + M 

2 . 0 2 +M^ 0 + 0 +M 

3. N0+M^N+0+M 


(20) 

(21) 

(22) 


and the bimolecular exchange reactions 
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+ I i NO + 0 

(23) 

+ 0 - WO + N 

( 24 ) 


where M denotes any one of the constituents that is catalytic in its effect 
on a reaction. 

The bimolecular exchange reactions (23) and ( 24 ) and the NO 
dissociation-recombination reaction ( 22 ) are considered to be the principal 
mechanisms that establish the balance of nitric oxide in the system. From 
equilibrium studies it is expected that the nitric oxide will be present in 
small enough amounts to have negligible effect on the compressibility Z and 
enthalpy h. However, its presence cannot be ignored since the foregoing 
reactions taken as a chain can significantly affect the rate of formation of 
N 2 from N atoms (ref. 11 ) . Since the amount of NO would be relatively con- 
stant in such a chain reaction, the steady-state approximation (in which the 
net rate of production of NO is taken to be zero) will be employed. This 
approximation leads to an expression for the amount of NO present in terms of 
the amounts of the more abundant species. That expression can then be used in 
the rate equations so that the amount of NO no longer appears as a variable. 

Upon completion of the foregoing steps and rearrangement, the rate equa- 
tions can be written 


<*(w Ne - w N ) 

( w W e - ( w 0 e 

1 

4 

O 

hT 

0 

+ dw Ng 

( 25 ) 

dt 

t kw 

t no 

dt 

d(w Qe - W 0 ) 

II 

teS 

0 

1 

c 

HD* 

fe! 

1 

0 

CD 

1 

w Q )fo 
w + 

dw 0 e 

(26) 


dt t on t 00 at 


Additional details of the derivation are given in the appendix. The quanti 
ties %, fo, and r are given in terms of the state variables by the 
relations 2 

(w N - 1)(1 - V N ) + (y N - 1) 

% = 




( 27 ) 


(wo e - l)(l - wo) + (yo - l ) 2 

fo i 

(w 0e - y 0 ) 

_i_ = _L 

t en t n ( T m) s 


(28) 

(29) 


1 

T oo 

1 

t no 

1 

t on 


1 _ 

T o 


+ 


1 

( T od 


“T k f 4 k f 5 ^ 
§ kf 4 kf 5 (l - 




w 


0 , 


) 


) / y N - y N \ 2 

\ v N e - %/ 
/wp - y 0 y 

vo e - yoj 


(30) 

(31) 

(32) 
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_1_ 

t n 


(l + a)% 


kf l { 1 2 ^ + " y N) 6 l + ^N 5 ! 


+ & 


W 0 o°2 ^ , v R 0 , -0+ 

5 X + (wq - Yq)^! + Yo^i 


( 33 ) 


T 0 (1 + a)% f sl 2 


.{■ 


0 , J 1 ■ % , ' Xc N c n + 

— — — ■{ £L S 2 + (wjf - y w )s 2 + y N 5 a 


+ ^ 1 g — + (w 0 - y 0 ) & 2 + yo s 2 


} 


(3*0 




kf 3 p 


[Kq(1 + a)% 


1+ Wjj + ^ (l+ Wq) 


(w 0 - yo) 


(35) 


1 ^ 

( T nn)s G 


•^fsd - w N e )(J 


W N " 


kf 3 P 


A\ - W L%(l +a )%J L N b 


1 + v w + ^ (! + 


»o> (»w - y w )} 


( 36 ) 


0 - (1 + sQmJ 1 ' W ° a>(v ° : yp) k fi + (1 ~ V (v * : yM) k, 
L (» 0e - y 0 ) a 4 (^ e - y H ) a 

kf +^a)ii]| 1 + " I * + I* 1 + ' w o ) ] } 


.Ks(l 


(37) 


Kg = 3-98! 


•500 


0.15 


eX P 


75490' 


(38) 


The equilibrium constant, Kq, was obtained from reference 8 . The remaining 
undefined functions and constants are given in table II. The rate constants 
utilized here are taken from reference 12 . 

The quantities and to are the values of and tqo> respectively, 

that would apply in the absence of the reactions involving NO (i.e., for pure 
nitrogen or pure oxygen) • For those cases, the coupling terms in equa- 
tions ( 25 ) and ( 26 ) would be absent. Thus the quantities t 0N> ( t Nn)s> 

and ( T Oo)s can " be associated with the reactions involving NO, while and Tq 
are due to the nitrogen and oxygen dissociation reactions. It is Instructive to 
observe the behavior of these quantities evaluated under equilibrium conditions. 

Figure l(a) is a plot of the quantity pT^ at equilibrium as a function 
of temperature with density as a parameter. The reaction times for dissocia- 
tion of pure nitrogen and dissociation of pure oxygen tq also appear in 


9 



the figure. The products pT are plotted rather than t because of the 
resulting suppression of density dependence. At low temperatures is 

approximately equal to Tq, but with increasing temperature it begins to 
exceed t q and eventually approaches Tjj. Thus in the low temperature range 
the nitric oxide chain reactions accelerate the dissociation (or recombination) 
of nitrogen in such a manner as to make tjjjj more nearly equal to tq than 
to T|j# In addition equations (25) and (26) indicate that the NO reactions 
bring about a direct coupling between the degrees of dissociation of nitrogen 
and oxygen such that one cannot be in equilibrium unless the other is also. 

The characteristic time t^q for the rate of dissociation of nitrogen due to 
a nonequilibrium condition of the oxygen is plotted in figure l(b) . 

Figure 2(a) is a plot of ptqq. It is seen that tqq does not deviate 
greatly from the reaction time of pure oxygen tq, particularly at the lower 
densities. However, equation (26) shows that a nonequilibrium condition of 
the nitrogen will affect the rate of dissociation of oxygen. The reaction 
time' corresponding to this effect tq^ is plotted in figure 2 (b). Note that 
the values of the reaction times appearing in figures 1 and 2 were evaluated 
at equilibrium. However, their dependence on the composition is weak and will 
not deviate greatly in nonequilibrium states. This is not true of the factors 
f n and ±q appearing in equations (25) and (26) and evaluated in equations (27) 
and (28). These quantities are usually of order 1.0 but can become large in 
some circumstances . For example, equation (27) shows that when the local 
equilibrium value of the degree of dissociation of nitrogen (w^ - yjj) becomes 

small, the factor fu will be large. For this reason, if the degree of dis- 
sociation of nitrogen is small, the nitrogen reaction can be nearly in equi- 
librium while the oxygen reaction is frozen. This is an exception to the 
foregoing argument to the contrary that was based on the values of t without 
consideration of the values of f]\[ and f q. 


NONEQUILIBRIUM CHANNEL FLOW CALCULATIONS 


The relations presented in the previous section are applicable to a wide 
range of flow problems including steady and unsteady one-, two-, and three- 
dimensional flows. We shall consider only the quasi -one -dimensional case here. 
In previous calculations of this type one difficulty that has plagued inves- 
tigators is the singular perturbation effect that arises in near equilibrium 
flow. The nature of the difficulty in the case of one-dimensional flow can be 
clarified by reference, for example, to equation (l6). In a numerical inte- 
gration procedure such as the Runge-Kutta method, the right side of equa- 
tion (l6) would be evaluated numerically several times at each step. In a 
near equilibrium region the factor (ty^)" 1 becomes large and e Y ± = ^ve± ~ e vi 
becomes small while the product retains an intermediate value that cannot be 
neglected. From a numerical viewpoint it is clear that evaluation of the 
product will entail a small difference when £ve± “ e vi is small. One 
observes that the chemical rate equations (eqs.(25) and (26)) contain the same 
near equilibrium indeterminancy. For a gas system involving many reactions 
the problem is amplified by the occurrence of many terms of this type which 
may exhibit improper behavior at different points in the flow field. Various 
remedies for the difficulty have been found (refs. 26 to 28, 11, 19, and 20). 
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It is perhaps fair to state that in practice these methods lead to a very- 
small step size in the numerical integration. In contrast, the procedure to 
be described here leads to accurate integrations that require no smaller step 
size than is needed for the corresponding equilibrium flow. Furthermore, the 
same nonequilibrium equations are used in the entire flow field including non- 
equilibrium, near-equilibrium, and equilibrium regions. Recently three new 
methods for integrating chemical rate equations have appeared in the litera- 
ture (refs. 29 to 3l). These methods also permit greatly increased step size 
in near- equilibrium flows. We have developed the required technique for inte- 
grating an arbitrary number of equations (ref. 32 ), including coupled- 
vibrational nonequilibrium, but the procedure for the general case will not be 
considered here. 


The crux of the method lies in the integration of equations (16), ( 25 ), 
and (26) over a small interval At. For example, equation (16) may be inte- 
grated exactly to obtain 


e v± 


2 



pt 2 dr} 

Ai T vi 


A V i T vi exp(u)du 



( 39 ) 


The subscripts 1 and 2 define the interval of integration and At = t 2 - tx* 
Using mean value integral theorems, we may write equation (39) 


G vi = e vi exp| 

2 1 


At 


( Tv i)i x . 


+ (AviTvi) g 2 -jl - exp[- 


(40) 


The subscripts |x and denote that Tvi and the product (AviTvi) are 
evaluated at different instants of time in the interval At which includes 
both £ x and | 2 • The closed form integration given by equation (40) is exact, 
but in practice it is difficult to evaluate since the quantities (Tvi)^ and 
(Avi T vi)| 2 are not readily obtained. However, if T V i and A v q are evaluated 
at the midpoint of the interval At and are monotonic over the interval, it 
can be shown that the error in e V i is of order 


evi(At ) 3 exp ^ (4l) 

In nonequilibrium regions of a channel flow calculation the Tvq has 
small enough values that the exponentials in equation (40) could be expanded 
to order (At) 2 without loss of accuracy. There would then be no difficulty 
in a numerical integration of equation (l6) or for that matter of equa- 
tions ( 25 ) and (26). For near-equilibrium regions, however, t v i becomes 
small, so the expansion is valid only for very small At. This is the reason 
a small step size would be required for numerical integration of equation (l6). 
If the closed form integration is used instead, the error becomes vanishingly 
small under these conditions because of the exponential factor in the error 
term above. 
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The ideas incorporated in the integration of the vibrational rate equation 
are adaptable to integration of the coupled chemical rate equations given by- 
equations (25) and (26). For that purpose these equations can be written in 
the form 

dei . , . . 

-jrr ~ -an€x + a i2 e 2 + Ai (42) 



where 

ei = 
€ 2 = 
an = 
a 12 = 

a 2 i = 

a 22 = 

Ax = 
A 2 = 


L — a 22 £ 2 + A 2 

(^3) 

w N e " W N 

(44a) 

O 

CD 

1 

O 

(44b) 

f N/ T M 

(44c) 

£ 'o/ t no 

(44d) 

f w/ T ou 

(44e) 

f oA 00 

(44f ) 

^- w N e 

(A4g) 

# 

a> 

(44h) 


The quantities €3. and € 2 are chemical nonequilibrium variables similar in 
character to e V j_ in that they are measures of the departure from equilibrium. 
The closed-form solution of these equations is 


6i = Bii exp(-AiAt) + B i2 exp(-A 2 At) + Cq (i = 1,2) (45) 


where 



7|"j( a ii + a 22 ) ± [(an - a 22 ) + 4ai 2 a 2 i] 


i/; 


B n = [(an - A 2 )€p i - ai 2 €p 2 ]/(Ai - A 2 ) 
B12 = -[(an - Ai) 6 p^ - ai 2 €p 2 ]/(Ai - A 2 ) 
Ci = (a 22 Ai + ai 2 A 2 ) /( ana 22 - ai 2 a 2 i) 
e P x = ( £ i)o “ 


(46a) 

(46b) 
(46c) 
(46d) 
( 46e) 


The quantities B 22 , B 2 i, C 2 , and ep can be obtained from the above relations 
by interchanging indices 1 and 2. The quantities (ei)o and (e 2 )o are the 
values of £i and e 2 at the beginning of the interval At. This is not an 
exact integration since the aqj and Aq are evaluated at the midpoint of the 
interval At in the same manner as suggested in the discussion following the 
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integration of the vibrational rate equations. The error is of the same order 
as given by the relation, equation ( 4 l). 

The remaining equations for quasi-one-dimensional flow can be found in 
many of the foregoing references (e.g., ref. 8). To save space these will not 
be repeated here. The closed form integration of the rate equations could be 
combined with the remaining flow equations in several ways. We have used an 
iteration procedure. Our program includes options for performing integrations 
with a specified pressure distribution p(x) or a specified area distribution 
A(x) . In the latter case a procedure is included for iterating to find the 
mass flow in the presence of nonequilibrium effects ahead of the throat. Inte- 
grations from the stagnation region to the throat are then usually repeated 
three or four times. Provision is made for iterating to find the conditions 
behind a normal shock, and by setting A(x) equal to a constant, the nonequi- 
librium flow behind a normal shock can be calculated. 

Our program contains options for including or excluding uncoupled vibra- 
tional nonequilibrium effects. Any of the rate constants can be separately set 
equal to a large value to obtain calculations corresponding to partial or com- 
plete equilibrium. Also the rate constants can be set equal to zero to obtain 
frozen flow calculations. In all of these cases the integration method is 
unchanged and remains the same in all regions of the flow since it is valid in 
the entire spectrum of equilibrium, none qui lib rium, and frozen flows. The step 
size is variable to conform with the requirement of a specified permissible 
deviation of the values of the key variables computed in two different ways 
(direct calculation and extrapolation from results at previous points) . Thus 
nonanalytic behavior of the specified function A(x) leads to a decrease in 
step size and is followed faithfully by the calculation. Accurate calculation 
of a nozzle flow (less than 0.1 percent change from reducing the step size 
criterion by a factor of 2 . 0 ) typically requires less than 200 steps. The 
average computing time on an IBM 7 ^ 9 ^ is about 45 seconds per case for equilib- 
rium, nonequilibrium, or frozen-flow calculations. 


EVALUATION OF APPROXIMATE MODEL 


The major consideration here is to evaluate the use of the steady-state 
approximation as a means of simply including the effects of the bimolecular 
exchange reactions. This is perhaps best done by comparison of the results 
with other numerical calculations and with experiment. For the comparisons 
with other numerical calculations we have chosen the work of Emanuel and 
Vincenti (ref. 9) and the work of Eschenroeder, Boyer, and Hall (ref. 11 ) . 
Complete listings for the latter calculation including temperature, pressure, 
and composition profiles were kindly provided upon request by Dr. D. W. Boyer 
of the Cornell Aeronautical Laboratory. The experimental comparisons are 
based on references 21 to 23, 33 j and 3^* 

Figure 3 includes temperature profile plots of a nozzle calculation 
obtained from the model and from Emanuel and Vincenti (ref. 9). The small 
differences indicate: (l) that there are no gross errors in either calcula- 

tion; and (2) that the steady- state approximation for the nitric oxide chain 
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reactions is valid for this case. The other flow field variables such as pres- 
sure, density, and velocity are less sensitive to the gas model used than is 
the temperature* No observable differences were noticed between the published 
curves and the model results for these variables. 

The compositions resulting from the two calculations are compared in 
figure 4 . Here the differences are greater because of the omission of nitric 
oxide from our model. Consequently, the amounts of N 2 , 0 2 , N, and 0 are all 
larger in our calculation than in that of Vincenti and Emanuel. However, the 
situation here is similar to that which exists in the comparisons at equilib- 
rium, previously discussed. Namely, replacing NO with varying amounts of the 
other species in the equations of state has less effect on the enthalpy h and 
the compressibility Z than might at first be expected. This is borne out by 
the close agreement between temperature distributions shown in figure 3 . The 
composition profiles may be computed more accurately if a correction discussed 
in the appendix is used. This was not done here to emphasize the fact that a 
good comparison is obtainable without accurate values for the concentration 
variables . 

Figure 5 shows four separate plots of temperature versus area ratio. One 
pair of curves represents a comparison of the model with the results of 
Eschenroeder, et al. (ref. ll) where the chemistry is in nonequilibrium and the 
vibrational degrees of freedom are in equilibrium. The other pair of curves is 
a comparison of calculations assuming complete equilibrium. One observes that 
the equilibrium calculations are in excellent agreement, as would be expected. 
The differences are greatest at the lower temperatures (and far downstream of 
the throat) where the omission of nitric oxide from the equations of state has 
its largest effect on the enthalpy and compressibility. The differences are 
greater, however, in the comparisons of the nonequilibrium calculations 
although the agreement is still excellent at the throat where the flow field 
is in local equilibrium. Downstream the flow field remains in equilibrium 
longer for the model than the results of Eschenroeder, Boyer, and Hall indicate. 
Whether the differences are due to the NO approximation or to the different 
methods used in integrating the near-equilibrium region of the flow field is 
difficult to determine. We did compare the nitric oxide concentration obtained 
from the steady- state approximation with that obtained from a more accurate 
relation (see appendix) and noted little difference. 

In figure 6 results are compared for the corresponding pressure profiles. 
Although the differences are smaller they still show the same character as 
discussed above, that is, the comparisons are in good agreement in the equilib- 
rium region of the flow field near the throat. Figure 7 shows a comparison of 
the composition profiles. Here the concentration variables are corrected by 
the method described in the appendix. The curves corresponding to the model 
are labeled 7i e or y± depending on whether the variables result from an equi- 
librium or nonequilibrium calculation. The results of Eschenroeder are indi- 
cated by ( ) c . As before, the agreement is excellent near the throat and 

deviations begin in the near-equilibrium region of the flow field. 

Nagamatsu and Sheer have published experimental data on vibrational relax- 
ation and recombination of nitrogen and air in hypersonic nozzle flows 
(refs. 33 and 34 ). Their data are given in the form of the ratio of static 
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pressure to reservoir pressure versus reservoir temperature at various speci- 
fied reservoir pressures, and at a station in the nozzle where the area ratio 
(a/a*) is equal to 144. We have made a series of (about 180 ) nonequilibrium 
channel flow calculations for the nozzle shape specified in reference 33 . The 
results are shown in figures 8 and 9 for nitrogen and air, respectively* In 
each of the figures the circles are the data points of Nagamatsu and Sheer 
while the four solid curves represent calculations using the model. For brev- 
ity the solid line curves are labeled with the numerals (l) through (4), respec- 
tively, designating the following types of flow field calculations: 

(1) Complete equilibrium 

( 2 ) Chemical nonequilibrium; vibrational equilibrium 

( 3 ) Chemical and vibration nonequilibrium (for pure nitrogen (figs. 8 (a), 
(b), and (c)), the vibrational rate constants are from ref. 19 , while 
for air (figs. 9 (a), (b), and (c)), the vibrational rates are those 

of ref. 35 ) 

(4) Frozen (at the reservoir conditions) 

The equilibrium and frozen flow calculations are included to establish the 
trends Indicated by the experimental results. Equilibrium flow field calcula- 
tions correspond to an infinitesimally small relaxation time (infinite rate 
constant), while frozen calculations correspond to an infinitely large relaxa- 
tion time (zero rate constant). 

The model results for pure nitrogen are given in figures 8(a), (b), and 
(c) and represent exact calculations since no assumptions are involved (the 
steady-state approximation does not apply here). One observes that below 
5000 ° K there is a negligible amount of nitrogen molecular dissociation since 
curves (l) and ( 2 ) are coincident. The nonequilibrium effects indicated by 
( 3 ) and the experimental data are then due primarily to vibrational nonequilib- 
rium. The experimental data points in figures 8(a), (b), and (c) all lie 
between the curves labeled ( 2 ) and ( 3 ) indicating that the vibrational rate 
constants used here are too small. This is in accord with references 21 and 22 
where It is suggested that the vibrational rates obtained from shock tube 
experiments must be greater by a factor of 15 in order to agree with expanding 
flow measurements. 

In figures 9(a), (b), and (c) the experimental results are compared with 
the model calculations for air. The experimental results are, in general, on 
or above ( 3 ), the curve representing vibrational equilibrium. The trend of 
these comparisons for air is in accord with references 22 and 23 where it is 
shown that the chemical equations given by equations (20) to (24) with vibra- 
tions in equilibrium (i.e., the model given in ref. 11) provide the mechanism 
for experimental agreement. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., May 5; 1965 
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APPENDIX 


SUPPLEMENTAPY EQUATIONS AND DERIVATIONS 


Relationship Between Lighth.il! Variables and Species Concentrations 


The species concentrations (moles of i per unit mass of gas) of the 

nitrogen components are given by the relations 


7 n 2 



(Al) 


7N “ ( W N ” y N )7t,N (A 2 ) 

7 N + = YN 7 t,N ( A 3) 

Corresponding expressions for 7q 2 ? 7q y and 7 q+ can be obtained from the above 
relations by replacing the subscript N with an 0, The values of 7 e - and 
the constants 7t,N> 7t,0 are 


7e~ = + 7 0 + = UN 


+ ? y ° 




(A4) 


7t,N = 


1 

(l + a)%[ 


(A5) 


7 t,0 “ § 7 t,N 

where ]% is the molecular weight of atomic nitrogen. 


(A 6) 


Chemical Rate Equations for Dissociation Reactions 


The rate equation for reaction 1 (eq. (20)) can be written (see, e.g., 
ref. 12) 

J / 

7 ifc " k“ 7 N 2 

The subscript 1 in dx/^/dt indicates that this is the rate of change of 
due to reaction 1. The index l refers to the catalytic body involved 
(e.g., N 2 , 0 2 , etc.). Upon substitution of equations (Al) and (A2) , this 
becomes 



(AT) 


dw-, 

dt” 



w, 


N 


7 t,N^ w N " 


(A8) 
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An expression for the equilibrium constant Kx in terms of the local equilib 
rium variables can be found upon noting that the quantity in square brackets 
is zero at equilibrium; namely, 

2p = 1 ~ v Ne 

Ki 7t » H " - y H ) 2 

Substituting this into (A8) and rearranging yields 

div N _ (]% ~ 
dt t n 

where % and t n are given in equations ( 27 ) and ( 33 )* The corresponding 
rearrangement of the rate equation for reaction 2 is formally the same, but 
with the subscript 1 replaced by 2, and the N subscripts replaced by 0. The 
l 

rate constants k-^ are given by 

kf ± = kf.S? (All) 

where the kf i and are listed in table II, which is based on reference 12. 


(A9) 


(A10) 


Steady-State Approximation 

Rate equations corresponding to reactions 3> and 5 can be expressed as 


-dgjj = ^0 = ilZNO = pk f / 7 0 2 7N - k Wo) 


dt dt dt 

d5y N _ _ d 57p __ d 5 7 N0 
dt 


o' (J u-b/NQ . / 1 \ 

— "if = f%s(7N a 70 - ^ 7N0?H) 


dt 


K 5 


(A12) 

(A13) 

(A14) 


/ 


The equilibrium constant Ka is given in equation ( 38 ) , while K 4 and K 5 can be 
expressed in terms of Kx, K 2 , and K 3 by the relations 


K _ K 2 
4 


K 5 = ^ 


Ka 


(A15) 

(A16) 


From the above expressions, the total rate of change of 7 N q can be written 
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(A17) 


d-^NO 

dt 



k f + 


k 4 


+ 



7 N0 




+ k f 4 7 0 2 7 N + k fs 7 N 2 7 0 


In the classical steady-state approximation (ref. 36), d7j\jo/<3.t is set equal to 
zero so that a value of 7^0 can be obtained from the above relation; namely, 




+ k f4 7 0 2 % 


+ k f5% 2 7 0 


k f 3 7 Z 


+ + k f 5 7p 

K 4 k 5 


(A18) 


where the superscript zero denotes that 7 jjq is a zeroth order approximation 
for the amount of nitric oxide. This value of 7^ 0 is then substituted in 
equations (A12) to (Al4) to obtain approximate rates of change of 7^ and Jq 
due to the reactions involving NO. Upon further substitution of equations (A.15), 
(Al6), (Al) , (A2), (A9) , and the corresponding relations for Jq , 7p, K 2 , and 
rearrangement, the rate of change of wp- due to the NO reactions can be 
written 


d.5% _ (wn s ~ vjj)fu (wp e - Wp)fp 

dt dt dt “ ( t NN) s t N0 


(A19) 


In the derivation of this equation, terms of the type d 4 yjj/dt appear. It can 
be seen that such terms are zero by noting for example that the foregoing term 
is proportional to d 4 7j^+/dt . The latter quantity is clearly zero (i.e., the 
rate of change of 7 jj+ due to reaction 4 is zero) . Expressions for the func- 
tions fp, (TjjjjJg, Tjjp are given in the text. Finally, equation (25) is 

arrived at by adding diw^/dt to obtain the total rate of change of w^. 

The corresponding relations for dwp/dt can be obtained from the above by 
interchanging subscripts 4, 5, and N, 0. 


Correction for Species Concentration Resulting 
From Nitric Oxide Omission 

Equation (Al8) yields the nitric oxide concentration resulting from the 
use of the steady-state approximation. To assess the accuracy of this approx- 
imation, 7^ may be compared with 7 |Jq obtained by integrating the nitric 
oxide rate equation given by equation (AI7) . This is done by using a method 
similar to that used to obtain equations (40) and (45). Equation (A17) may be 
represented by 

¥2° = -A 7h0 + B (A20) 
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The variable B includes all the terms in equation (A 17 ) that do not include 
7 ^q and the variable A represents all the terms with 7^0 factored out. 
Equation (A 20 ) is then integrated by the same procedure discussed in the text 
to obtain 

= (^no)o ex P(” A At ) + [1 " ex P( A At )] 7 ^o ( A 2 i) 

Here the ratio B/A is replaced by 7 ^q obtained from the steady-state approxi- 
mation and (7^0)0 represents the value of 7^0 at the beginning of the inter- 
val (see paragraph in text preceding eq. ( 4 l)). 

The values of A and 7^0 are obtained by evaluation at the beginning of 
the interval. The contained in A and 7 §q are corrected in the following 

manner. The corrections to the number of nitrogen and oxygen atoms, 7^ ^ and 
7t,0> are given by 

7 t,N = 7 t,N + 7 N 0 {& 22 .) 

7 t ,0 = 7 t ,0 + 7 N 0 ( A23 ^ 

where 7^ ^ and 7t 0 are given by equations (A5) and (A6), respectively. Sub- 
stituting 7 t,N into equations (Al) to (A3) and q into the corresponding 
relations for the oxygen species yields the corrected values of the species 
concentrations . 
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TABLE I.- FIXED CONSTANTS REQUIRED FOR NITROGEN AND OXYGEN MIXTURES 
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TABLE II.- RATE CONSTANTS 


(a) Chemical N one qui lib r ium Rate Constants 


Reaction 

Catalytic 

Forward rate 

Efficiency 

equation 

body. 

constant 

factor. 

number, i 

l 

k f,i 


1 

n 2 

3.0X10 21 T" 3/2 exp (hm260^ 

1 


0 


l/3 


02 


1/3 


H 


5 

2 

O2 

3.6xi0 21 T" 3/2 exp 

1 


0 


0.59X10' 3 T 


n 2 


1/3 


N 


1/3 

3 


! 5 • l8xl0 21 T" 3/ 2 exp(tS^2i) 


4 


1.0X10 12 T 1/2 exp (^^ 2 °j) 


5 


4.5X10 13 exp(^8|li) 



^-Depending on reaction, kf^i has dimensions of cm 6 /mole 2 -sec 
or cm 3 /mole-sec . 


(h) Vibrational None equilibrium Rate Constants 

pTy = 1 . 6188x10 " 3 exp( 101. 44 /t 1/3 ) dynes -sec/cm 2 

pry^ = 1.11531X10 5 nTt exp(l5^ .0/T 1/3 )dynes-sec/cm 2 
where p is the local value (or static) pressure 



Temperature, T, 1000° K 


(a) PTjjip Ptq, pipj versus temperature. 

Figure 1.- Product of density and the characteristic relaxation times that occur in the rate equation 

for nitrogen. 



1 o ^IO(/ 5T NO) 


-8 


-10 


-12 



! 


-14 


0 


5 10 

Temperature, T, 1000° K 

(b) pt^q versus temperature 


Figure 1.- Concluded. 



Temperature, T, 1000° K 

(a) PTqq, pT Q , pt k versus temperature. 

Figure 2.- Product of density and the characteristic relaxation times that occur in the rate equation 

for oxygen. 
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(b) pt on versus temperature. 


Figure 2.- Concluded. 




Composition, y° , moles per gram of gas 


I O ' 1 - Air 

n p r = 1248 atm 

T r = 7950° K 




Temperature, T, 1000 



Figure 5.- Comparison of equilibrium and nonequilibrium nozzle temperature profiles obtained from 

Eschenroeder et al.. with those of the model. 



Figure 6.- Comparison of equilibrium and nonequilibrium nozzle pressure profiles 
obtained from Eschenroeder et al.. and the model. 
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Reservoir temperature, T r , 1000° K 

(a) Reservoir pressure equal to 6.8 atm (100 psia) . 

8.- Comparison of experimental results for pure nitrogen obtained from Hag ama tsu and Sheer 
(ref. 33) and from a calculation using the integration techniques of this report. 
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Figure 8.- Concluded. 
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(a) Reservoir pressure equal to 6.8 atm (100 psia) . 

Figure 9 .- Comparison of experimental results for air obtained from Nagamatsu and Sheer (ref. 33 ) with 

the air model. 
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(b) Reservoir pressure equal to 13.6 atm (200 psia) . 
Figure 9 .- Continued. 
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